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Abstract 

In non-minimal Higgs mechanisms, one often needs to minimize highly symmetric 
Higgs potentials. Here we propose a geometric way of doing it, which, surprisingly, is 
often much more efficient than the usual method. By construction, it gives the global 
minimum for any set of free parameters of the potential, thus offering an intuitive un- 
derstanding of how they affect the vacuum expectation values. For illustration, we apply 
this method to the S4 and ^4-symmetric three-Higgs-doublet models. We find that at 
least three recent phenomenological analyses of the ^4-symmetric model used a local, not 
the global minimum. We discuss coexistence of minima of different types, and comment 
on the mathematical origin of geometrical CP- violation and on a new symmetry linking 
different minima. 

1 Introduction 

Finding out the nature of the electroweak symmetry breaking (EWSB) is one of the hottest 
topics in high-energy physics these days. The first LHC data on the Higgs-Uke resonance at 
126 GeV show intriguing deviations from the Standard Model (SM) expectations, [T], and 
many believe that they hint at a non-minimal Higgs mechanism of EWSB. 

In the past decades, many non- minimal Higgs sectors have been considered, [2]. Typically, 
these sectors involve several Higgs fields interacting via the scalar potential, which is often 
invariant under a group of Higgs-family transformations. Once the potential is written, one 
then proceeds by minimizing the potential, finding the vacuum expectation value (vev) align- 
ment, expanding the potential near this point, and calculating phenomenologically relevant 
quantities. 

The standard procedure for minimization of the potential is to parametrize the vev's via 
(possibly complex) Vi, calculate V{vi), then set all dV/dvi = 0, solve these equations for Vi, 
and finally check that the hessian at this point is positive definite. Alternatively, one can start 
with the desired vacuum configuration and build the potential with a prescribed symmetry 
around the vacuum point. 

Although this method usually works well, there are several reasons why one might not 
be completely satisfied with it. First, sometimes the equations cannot be solved analytically. 



Second, if they are solvable, they do not always give a clear intuitive picture of how the vev 
alignment depends on the free parameters of the potential. Third, the potentials can support 
several local minima, and in order to be sure that one works with the global minimum, one 
must check that all other possible minima of the given potential lie above the chosen one. 
Unfortunately, this check is done very rarely, merely because this is a difficult task on its own. 

In this paper we propose a geometric method of minimization of potentials which is free 
from these drawbacks. Its advantages are: 

• by construction, it gives the global minimum of the potential, 

• one does not need to repeat the calculations for various regions of free parameters; once 
the main geometrical object (the orbit space) is constructed, one gets answers for all 
allowed free parameters, 

• even in cases when the analytic minimization is impossible, it can still give useful infor- 
mation e.g. number of degenerate minima, symmetry breaking properties, etc. 

The main drawback of this method is that it is handy only for highly symmetric potentials, 
so that the number of free parameters of the potential is small. However this situation takes 
place in many particular realizations of non-minimal Higgs sectors, and this method should 
indeed be useful in practical calculations. 

In this paper, we will illustrate this method with the particular versions of the three- 
Higgs-doublet model (3HDM), the ones with the and S4 Higgs- family symmetry group. 
This models has been actively studied in the past few years [21 IH El El E] , with the idea that 
it might provide a natural explanation to the patterns observed in the fermion mass matrices. 
Minimization of the potential was conducted in these papers in the standard way, and several 
vev alignments were used. However, we will show below that some of them correspond in fact 
to a local, not a global minimum; a fact that apparently went unnoticed up to now. 

The structure of the paper is as follows. In Section [2] we describe our main idea, which 
we formulate for convenience in the context of multi-Higgs-doublet models. We illustrate the 
general construction in Section [3} where we discuss in detail the 5*4 and A4-symmetric 3HDM. 
In Section |4] we comment on previous publications, discuss various additional aspects of the 
method, and finally draw our conclusions. In Appendix, for completeness, we present the 
Higgs mass spectra for all global minima possible in 5*4 and A4-symmetric 3HDM. 

2 Geometric minimization of symmetric potentials 

Although the method we propose is rather general and can be applied to a broad range of 
extended Higgs sectors and perhaps beyond, we prefer to expose it in the context of iV-Higgs- 
doublet models (NHDM). This will allow us to keep the notation simple and, at the same 
time, get prepared for the particular applications in 3HDM. 

2.1 Orbit space in NHDM 

In the iV-Higgs-doublet model we introduce Higgs doublets (pa, d = 1, ■ ■ ■ ,N with identical 
electroweak quantum numbers. The general renormalizable Higgs potential of NHDM is con- 
structed from the gauge- invariant bilinear combinations (0^06) [lOj, which describe the gauge 
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orbits in the Higgs space 15 The space of gauge orbits (the orbit space) can be represented as 
a certain algebraic manifold in the space of these bilinears. It is convenient to group them in 
the following N"^ real bilinears: 



IN 



2N 

a a,b 

where A* are the generators of SU{N). The orbit space as a manifold in the euclidean space 
of these bilinears was characterized algebraically and geometrically in [8] . It lies between 
two forward cones defined by 

^o>0, ^^^~_\^ rl < < rl (2) 

In addition, neutral vacua always lie on the surface of the outer cone = Tq, while charge- 
breaking vacua lie strictly inside, -r^ < Tq. 

In the formalism of bilinears, the Higgs potential takes the form of a general quadratic 
form of ro and rf. 

1 1 
V = -Moro - Min + -Aooro + Ao^rori + -A^jrirj . (3) 

The minimization of the potential can then be cast into a geometric condition of contact of two 
algebraic surfaces constructed in : the orbit space and the equipotential surfaces defined 
hj V = const, see details in [9]. 

In the particular case of 3HDM, the bilinears are 



2 ' " 2V3 

ri = Re(0|02) , = Re(0j50i) , tq = Re(0^03) , 

ra = Im(0j02) , rg = lm(0^0i) , rj = Im(0^03) . (4) 
The orbit space in 3IIDM is defined by 

ro>0, <rl, y^dijkTirjrk = — -Vq , (5) 

and the modulus of the vector f is restricted as 

\rl <^< rl . (6) 



2.2 Geometric minimization: the main idea 

The crucial feature of passing from fields to bilinears is that the Higgs potential is simplified; 
it becomes a quadratic form of these new variables. This transition was used in the two-Higgs- 
doublet model (2HDM) [11], and it allowed one to observe and exploit interesting geometric 

"'^Strictly speaking, (j)a are operators acting on the Higgs Fock space; however, for the purposes of this paper 
one can view them as doublets of complex numbers. 
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features of the potential both in 2HDM \l2l and in multi-Higgs-doublet models [9]. Here we 
propose to go further in this direction, and introduce new variables in terms of which the 
potential becomes a linear function. 

To this end, consider an NHDM potential with a sufficiently high symmetry so that Mj = 
in ([3]). Absence of this term is a hallmark of so-called "frustrated symmetries" in NHDM which 
were discussed in p!3j. The quartic part of the potential contains k different terms, k usually 
being rather small for a highly symmetric potential. Let us generically write the potential as 

k 

Here Xi are the quartic terms divided by Tq, with xq = 1 by convention, and Aj are coefficients 
in front of them. 

Let us now consider the variables Xi, which can always be chosen real. Calculating them 
for all possible field configurations (or for all possible values of r's inside the orbit space) will 
fill a certain region in the space M.^. This region, which we denote by F, is the orbit space 
"squashed" into the Xj space. Note that the map from r's to Xi is not, generally speaking, 
injective because different r's can correspond to the same point x,. 




(a) (b) (c) (d) 



Figure 1: Two-dimensional illustration of the geometric minimization method. Shown are four 
cases with the same orbit space F (represented by the shaded region) but with different values 
of the parameters Aj, which define the direction of steepest descent n. The four cases show 
how the local geometry of F determines the stablity and degeneracy of the global minimum; 
(a): locally convex geometry leads to a single global minimum indicated by the dot, whose 
position is sensitive to the exact values of A^; (b): a minimum at a vertex is stable against 
variations of A^; (c): cusps separated by a concave region allow two distinct minima to coexist 
and be degenerate; (d): straight segments can lead to a continuum of global minima (shown 
by a thick line) for the special values of the parameters, implying presence of massless bosons. 



Suppose the geometric shape of F is known. Then the minimization of the potential 
proceeds in three simple steps. First, since the potential ([T]) is a linear function of Xj, we 
can introduce the "direction of steepest descent" of the potential, n = — (Ai ,. . . ,Ak). The 
potential can then be written as 

V = -Moro + rl (Aq - fix) . (8) 

Then the minimum of the potential is achieved at the points of F which protrude farthest in 
the direction of n. Once these points Xj are known, we can find their realizations in terms 
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of fields, and, finally, find the value of tq. Note that the positivity conditions require that 
Aq — nx > everywhere in F. 

Some phenomenologically relevant properties of the minima follow from this geometric 
picture, which we illustrate in Fig. [T] with a two-dimensional example. If F has a smooth 
and strictly convex local shape. Fig. [l^, then the minimum is unique in the Xi space. It 
is also seen that the minimum point continuously changes if parameters Aj are varied. If 
F, instead, has vertices, see Fig. then the minimum point becomes stable within certain 
regions of Aj variation (or alternatively, regions of possible directions of n). Note that such 
a feature is the origin of geometric CP- violation in multi-doublet models, see our discussion 
in Sect. 4.2 At the borders of these regions, two concurrent minima become degenerate 



and coexist. Fig. [T]:; crossing this border causes a first order phase transition between the 
two vacuum configurations. Finally, if F contains straight segments. Fig. [T]i, then for the 
borderline parameters Aj we get a continuum of global minima, which means that the model 
contains additional massless scalars. 

In this picture, the key object becomes the shape of F rather than the parameters of 
the potential. Once the symmetry group is fixed and F is constructed, many properties of 
the potential (points of minimum, their degeneracy and coexistence, patterns of symmetry 
breaking, the phase diagram of the model and phase transitions) can be immediately read 
from its shape. 



3 A4 and 5*4- symmetric 3HDM 



3.1 The potentials 

In this Section we will illustrate how the general method works with the example of A4 and 
5'4-symmetric 3HDM. 

The y44-symmetric 3HDM can be represented by the following potential 



V 



Ao 

3 



^'101 + 0202 + 0303 



A3 

+Ai 



(0l0l)^ + (0^02)' + (0503)' - (0l0l)(0502) - (0502)(0503) " (0503)(0l01 



+A2 
+A4 

or, in terms of bilinears 
V 



(Re0l02)' + (Re0503)' + {Re^l^iY 
(Im0l02)' + (Im0^03)' + (Im0],0i)^ 



(Re0{02)(Im0l02) + (Re0^03)(Im0^03) + (Re050i)(Im050, 



(9) 



-Moro + Aor^ + K^(rl + rl + rl) + K2{rl + + r^) + ^^{rl + r, 
+A4(rir2 + r^r^ + rgry) . 



(10) 



Here parameters Mq and Aj are assumed to take generic values. This potential is symmetric 
under the full achiral tetrahedral group Td isomorphic to x Z2 of order 24. This group is 
generated by independent sign flips of individual doublets, by cyclic permutations of the three 
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doublets, as well as by a specific type of generalized-CP transformation (the CP-conjugation 
combined with exchange of any two doublets). 

An alternative way to parametrize the potential was used in [H El [7] . Our coefficients are 
related with the coefficients of the alternative parametrization used in [7] as 

-^ = /i^, Ao = 3Ai + A3, A3 = -A3, Ai^2 = A4 ± As cose , A4 = -2A5sine. (11) 

If A4 = 0, we get the S'4-symmetric 3HDM. In the alternative parametrization, this is 
equivalent to setting e = 0. The potential becomes now symmetric under the full achiral 
octahedral group Oh isomorphic to ^4 x Z2 of order 48, which is generated by sign flips of the 
individual doublets, by their permutations, and by the CP- conjugation. 

Since the classiflcation of the finite realizable symmetry groups of the scalar sector in 
3HDM is now known [T3] , we know that restricting the parameters further will never produce 
any larger finite symmetry group. It can only lead to continuous symmetry groups, which 
are necessarily frustrated and must therefore be spontaneously broken and produce massless 
scalars [I3]. We disregard this situation on phenomenological grounds. 



3.2 The orbit space in the 6*4 case 

Let us start with the more restricted model, the S'4-symmetric 3HDM. Written in terms of 
bilinears, the potential takes form 

V = -Moro + Aorl + Ai{rl + rj + rl) + A2{rl + rl + r^^) + A^irl + rl) 

= -Moro + rl{Ao + Aix + A2y + A3Z) , (12) 

with the vector {x,y,z) playing the role of Xj. The positivity conditions for the potential 
require that 

Ao + Aix + Aay + Aa^ > (13) 

everywhere in the orbit space. Using the properties of bilinears mentioned in Section |2} 
we conclude that the three-dimensional orbit space F must lie inside the truncated pyramid 
defined by 

x,y,z>0, ^<x + y + z <1. (14) 
In addition, it turns out that y < 3/4. Indeed, y can be rewritten as 



(pi + P2+ p^y 



(15) 



where 

'>^^^>0, «:,, = p,p,-|0t0,|2>o, Q^^ = ^0-0^, (16) 



^a'ra) 



with a = 1, 2, 3 numbering the doublets and a, /3 = +, denoting the upper and lower compo- 
nents inside doublets. The largest value of y equal to 3/4 is attained when, first, all ^ab = 
which selects the neutral vacuum, and then when the lower components of the doublets sum 
up as = 0. 
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The exact shape of the orbit space which we found by numerical method^ is rather com- 
phcated, see Fig. |2| left. It has the form of a wedge with the edge ed. x + z = l,y = and the 
convex backside at small values of x + z. However if we focus only on the phenomenologically 
relevant case of neutral vacuum, then we limit ourselves only to one of its faces defined by 
X + y + z = 1. The rest of F corresponds to charge-breaking vacua and is disregarded. 



y 

D 


C 











B 




Figure 2: Left: sketch of the orbit space F of the S'4-symmetric 3HDM in the (x, ?/, z)-space. 
Right: the neutral orbit space in the (x, ?/)-plain. On each plot, the four dots A, B, C, and D 
mark the positions of the possible neutral global minima. 

This neutral part of the orbit space has the shape of a trapezoid, which can be established 
analytically using (14) and (15). We show it in Fig.[2| right, on the (x,y)-plane (with z defined 
as z = 1 — X — y). The vev alignments corresponding to the four vertices of the trapezoid are 
(here we give only the relative magnitude of the vev's) 

A: (1,0,0), 5: (1,1,1), C = (e*"/^ e-^"/^ ±1) , D = (e^"/^ e-^"/^ 0) ~ (1, i, 0) , 

(17) 

while the straight segments joining them are 

AB : (t;i,t;2,t;3) with all G M, BC : {e'^' , e'^\ e'^') , (18) 



AD : (f 1, iv2, 0) with all Vi G 



CD : (e^^e-^^r) with cos 2^ = 



Of course, in each case we allow for arbitrary permutations of the doublets. For example, 
vertex D corresponds to six degenerate minima (1, ±i, 0), (1, 0, ±i), (0, 1, ±i). 



3.3 Minimization of the 6*4- symmetric potential 



Applying the methods of Section |2.2[ we immediately conclude that the S'4-symmetric 3HDM 
can have only four types of neutral minima without producing massless scalars, which cor- 
respond to the vertices (17). Thus, we located all possible positions of the global minimum 



without the need to calculate any derivatives. 

It is also possible to obtain conditions on Aj which lead to a minimum at each of these 
four points just by looking at the orbit space. For example, the vev alignment of the type 
(1,0,0) becomes the global minimum, when A3 < and Ai, A2 > A3. When these conditions 
are satisfied, the point A indeed lies farthest along the direction n. In addition, the positivity 



condition (13) in this case implies that Aq + A3 > 0. 



^We generated one million points with random up and down components of the three doublets. For each 
point we calculated the values of x, y, and z, and then plotted all points. By looking at the resulting 3D 
scatter plot from different angles, we reconstructed the shape and drew the sketch in Fig. [2] 
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3.4 Unexpected symmetry of the orbit space 

The 5'4-symmetric 3HDM possesses a curious feature which could be noticed but would receive 
no explanation with the usual calculations. 

First, field content of the scalar sector, after the electroweak symmetry breaking, is the 
following: apart from the usual three "would-be" Goldstone bosons, we have two pairs of 
charge-conjugate Higges H^, and five neutral scalars. The oscillation mode in the direction 
of vev's will be denoted as h, while the other neutral Higgses are generically labeled as Hi. In 
(generalized) CP-conserving cases, these can be additionally classified as (generalized) CP- 
even and CP-odd states. 

Let us now calculate the masses of the physical Higgs bosons in the two vev alignments: 
(1, 1, 1) and (1, 0, 0). In both cases we use v"^ = vl + V2+ vj. The alignment (1, 1, 1) becomes 
the global minimum of the potential if 

Ai<0, Ao > |Ai| > -A2,-A3- (19) 
The minimum point is then parametrized as {v,v,v)/\/6 with 

P 1, . , 9 v^Mo lAil . , , , , 

ml.± = -\Ai\v — — - — j—— j- (double degenerate) , 

< 2 2 Aq — I A]^ I 

2 1/1 A I * N 2 V^Mq |Ai| + A2 , , , , , , 

iTT-Hi = 2^'^' ^ — 2 — Ao - [All degenerate) , 

1/1 A , A X 2 Mo lAil A3 , , , , , , 
-(|Ai| + Asjv = _ i^^i (double degenerate) , 

ml = l{Ao-\A^\y^^Mo. (21) 

The alignment (1, 0, 0) becomes the global minimum if 

A3<0, Ao > IA3I > -A2,-Ai. (22) 
Expanding the potential around the point (f , 0, 0)/\/2, we get 

2 1, . , 2 V^Mq IA3I , , , 1 , , . 

m„± — -\A3\v — — - — j— — j- (double degenerate) , 

'2 2 An — I Aa I 

mjj. = -(IA3I -|- A2)t'^ = " ^ '1 ^ ^ (double degenerate) , 

-(IA3I -|- Ai)v^ — ^ " ^ ^ (double degenerate) , 

ml = ^(Ao-|A3|)^2 = -^Mo. (24) 
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It is hard to miss a remarkable symmetry between these mass spectra: upon exchange Ai o A3 
they almost turn into one another. The only quantity that violates this otherwise perfect 
symmetry is the mass of one pair of neutral Higgses. 

The bizarre aspect of this almost perfect symmetry is that it is not a symmetry of the 
model. It would be a symmetry if there existed a transformation of fields that could swap 
X and z while keeping y unchanged. But such transformation does not exist. This is also 
consistent with the fact that Ai -H- A3 does not lead to an exact matching of the two Higgs 
spectra. 

We can trace the origin of this near symmetry from the shape of the orbit space F in the 
(x, y, z) space. Our numerical study offers very strong hints that this shape is indeed x -H- z 
symmetric; unfortunately, we do not have an analytic proof of this fact. Provided this is true, 
it explains why conditions ( [l9) ) and ( [22| and the charged Higgs masses (which are also related 
with the shape of the orbit space) are exactly symmetric. 

It is interesting to notice that, for both vev alignments, the spectrum of the Higgs bosons 
is 2HDM-like. Namely, we have only one value for the charged Higgs masses and three values 
for neutral Higgs masses, just as expected for the generic 2HDM. How this situations can 
be distinguished from the true 2HDM experimentally, and which observable quantities one 
should look at, is a separate issue worth investigating further. However the origin of this 
2HDM-like spectra is different in these two cases. In the vev alignment (1, 1, 1) it comes form 
the unbroken S'3-symmetry of the model, [3], while for the (1,0,0) alignment, the origin is 
the (9(2)-symmetry mixing the second and third doublets, which is manifest in the vev's and 
the mass terms. These two symmetry arguments are non-equivalent and cannot be related to 
each other. After all, the two vev alignments also differ in the number of degenerate vacuum 
points: four for the (1, 1, 1) and three for (1, 0, 0). 

The same relation holds between the other two possible vev alignments, see Appendix. 
It would be interesting to see if this approximate symmetry leads to other phenomenological 
similarities between these pairs of minima. 



3.5 The orbit space in the A4 case 



We write the generic potential of the yl4-symmetric 3HDM (10) in the way suggested in 
SecO 

V = -Mor-o + r^(Ao + Aix + A^y + Agz + A4t) , (25) 



with the same x,y,z as in Eq. (12) and with t = (rir2 + r^r^ + rQri)/rQ. 

Again, we focus on the neutral orbit space F, for which we choose x, y, and t as independent 
variables, and then z = 1 — x — y. The shape of F which arises from our numerical study 
is shown in Fig. [3] in the {x,y,t) space. Despite being rather complicated, it displays a 
remarkable triangle symmetry. Basically, it is a right circular cone oriented along the direction 
X — y = t = with the apex at the origin and with opening angle 7r/2 and with directrix lines 
of length 1. Parts of this cone starting from distance 3/4 from the apex are carved out. It 
has four flat faces which have the shape of deltoid (3- vertex cusped closed curve). These are 
located at the bottom face (deflned by z = 0) and at three side faces (one is y = 3/4 and 
the remaining two are obtained by the 27r/3 rotation of the cone). Directrix lines opposite to 
these three side faces have length 1 and extend up to the points B, C and C where three flat 
deltoid regions meet. The remaining portions of F are concave regions. The two-dimensional 
neutral orbit space for the 5'4-symmetric model. Fig. |2| right, is simply the (x, ?/)-projection of 
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Figure 3: Sketch of the neutral orbit space in the tetrahedral 3HDM viewed from two angles. 
Uniformly shaded regions correspond to flat faces, graded shading indicates curved faces. 

this F; correspondence between the vertex points in the two shapes should be clear. We only 
mention that point L>, which was a vertex in the 5*4 case, becomes a regular point on the rim 
of the cone in the A4 case. 

The conical shape of the orbit space can be understood in the following way. Let us 
introduce two real vectors 

a = — (ri,r4,r6) , ^ = — (r2, r5, ry) , (26) 
ro ro 

and denote the angle between them by ^. Then, x — cP,y — d^^t — (ab) = ^yxycos^. 

— * 

If it is possible for a given point {x, y) to find two parallel vectors a and 6, then the orbit 
space extends in the ^-direction up to i = y/xy, which precisely defines the cone. It turns 
out that parallel vectors a and h exist within the triangle < x + y < 3/4 (vev alignment 
is (f 1, f2e*", 0)), and along the straight segments defined by y = (vev alignment {vi,V2.iV^)) 
and defined by y/x = 3 (alignment (vie*'^/^, ^26"*'^/^, W3)). 

Let us also stress that the emergent triangle symmetry of the orbit space is not related to 
the symmetry of the potential, but is a feature of the orbit space itself. There is simply no field 
transformation that realizes rotations of the cone. In this aspect, this emergent symmetry is 
similar to the (x, z) reflection symmetry of the orbit space in the 5'4-symmetric model. 

3.6 Minimization of the yl4- symmetric potential 

The shape of the orbit space immediately leads to the list of possible phenomenologically 
acceptable global minima (i.e. minima not leading to additional massless scalars). These are: 
the apex of the cone (point A), the three vertices (points B, C and C"), and the rim of the 
cone (x + 1/ = 3/4, = xy). Any other point either leads to additional Goldstone bosons or 
is never a global minimum. 

In Appendix, we analyze all these points in some detail and give the Higgs mass spectra. 
Comparing these spectra for different minimum points also shows intriguing relation with the 
triangle symmetry of the orbit space. 
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To our knowledge, this is the first complete solution of the minimization problem in the 
yl4-symmetric 3HDM. 



4 Discussion and conclusions 

4.1 Comments on previous publications 

The yl4-symmetric 3HDM has received much attention, [31 HI El El [7j, because spontaneous 
breaking of A4 can generate interesting patterns in fermion mass matrices. In most papers, 
the authors just pick up a minimum with specific vev alignment, build a potential which 
indeed has a minimum at that point, and then proceed with analysis in the fermion sector. In 
particular, in IH E] the following vev alignment was used: 

(e^^e-^^r), (27) 

with r and a being independent real parameters. Fitting the model to the fermion observables 
gave r ~ 40 in and r ~ 240 in [5]. 

The authors of ^ aimed at a complete description of possible vev alignments in the minima 
of the 744-symmetric 3HDM. They critically reanalyzed the phenomenological situation arising 



at the minimum (27) and argued that it is strongly disfavored by the flavor physics constraints. 
However, they also considered this vev alignment as a viable solution of the minimization 
problem]^ 

This vev alignment is absent in our classification because it cannot be the global minimum. 
One can see it most easily precisely in the case of large r. Indeed, in our notation, this point 
corresponds to z being close to 1 and, therefore, it lies strictly inside the cone and close to its 
apex. This point would emerge at the surface of the cone only for angle a multiple of vr/S, 
but in this case one would get massless scalars. This means that even if this point happens to 
be a minimum, it is a local, not the global one. 

Thus, we prove that the phenomenological analyses of [HE] and, partly, of |7] correspond to 
a metastable electroweak vacuum. Apparently, this fact went unnoticed because the authors of 
these publications did not check whether the potential they got possesses a deeper minimum. 

It is generally believed that the Universe sits in the global minimum of the Higgs potential, 
because the early hot Universe while cooling down through electroweak temperatures either 
selected the global minimum or, even if it were stuck in a metastable vacuum, had enough 
time to tunnel to the absolute minimum. Taken seriously, this argument implies that much 
of the analyses of [H E] is not phenomenologically relevant on these grounds, even without 
appealing to the flavor physics constraints. 

Finally, let us note that the above remarks apply even to the softly broken symmetry, 
provided the symmetry breaking terms are not large. Indeed, by continuity arguments, the 
position and the depth of a minimum with the strictly positive hessian change continuously 
under a continuous variation of parameters. Therefore, for small y44-breaking terms the vev 
alignment (e*",e~*",r) still cannot be the global minimum. 



^ We are thankful to Luca Merlo who clarified to us the motivation behind the works . 
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4.2 Origin of geometric CP-violation 



The possibility for spontaneous CP- violation is one of the motivations behind studying multi- 
Higgs-doublet models. In this context it is often proposed not only that a Higgs-family sym- 
metry should allow for spontaneous CP-violation but also that it should stabilize the vev 
phases in the global minimum against variation of the free parameters. This situation is 
known as geometric CP violation, [151 [E] and was originally found in the A(27)-symmetric 
3HDM (though we note that the true Higgs-family symmetry group of that model is A(54)/Z3, 
see discussion in [T^). The relative vev phases arising in geometric CP- violation are called 
calculable because their values follow from group theoretic arguments and do not depend on 
the exact values of the parameters of the potentiaQ 

Using the method described in the present paper, we can pinpoint the mathematical origin 
of calculable phases in such models. They arise due to the presence of vertices in the orbit 
space r, see Fig. [T]D,c,d, or to be more specific, vertices at points corresponding to non-zero 
relative phases. Absence of geometric CP phases would imply convexity of the orbit space. 
Fig. [pi. So, it is not the symmetry of the model per se that allows for calculable phases but 
the choice of coordinates Xi selected by the symmetry, in which the orbit space has vertices. 

Our experience shows that the higher the finite symmetry group, the simpler is the ge- 
ometric shape of the orbit space F, and the more vertices linked by straight segments it 
possesses. This explains why it is natural that geometric CP violation starts to appear only 
for sufficiently large finite symmetry groups. 



4.3 Coexistence of different minima 

A priori, it might happen that, for some values of the parameters of the potential, two (or 
more) different types of the global minimum coexist and are degenerate. Fig. [1J3 illustrates 
this situation. Upon small variation of the parameters around this special point, one minimum 
point becomes the global minimum while the other turns into a local one, and it is clearly 
possible to make either of them the global minimum. This feature leads to a possibility of 
a first order phase transition upon smooth variation of the parameters, leading to important 
phenomenological consequences. It is therefore desirable to know, which models allow for such 
a possibility. 

It can be inferred from Fig. [T]^ that this can happen if the orbit space vertices separated 
by a concave region, that is, if it has cusps. If instead the orbit space is a convex body, this 
possibility is excluded. Possibility of a first order phase transition is, therefore, linked to the 
non-convexity of the orbit space. 

Our analysis shows that the orbit space of the S'4-symmetric 3IIDM is convex. Therefore, 
phenomenologically relevant global minima of different type cannot coexist in this case. 

In the y44-symmetric 3HDM, this possibility arises. Namely, generic points on the rim of 
the cone and one of the three vertices B, C, or C in Fig. [3] can be degenerate. It is also 
possible to make two among these three points degenerate, but not all three. Examples of 
such potentials can be readily costructed from geometric analysis of Fig. |3} 

However, we find the following statement from 116] inaccurate: "...t/ie calculable phase arising from geo- 
metrical CP violation is uniquely determined independently of the arbitrary parameters of the scalar potential." 
Indeed, parameters in any case must be sucli that the point realizing the geometric CP violation is the global 
minimum of the potential. This takes place only in a certain region but not everywhere in the space of the 
free parameters. 
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4.4 How general is the proposed method? 



In which cases does the geometric minimization method proposed in this paper become useful? 
Strictly speaking, it has no intrinsic limitation. For example, in the context of the multi-Higgs- 
doublet models one can start with an absolutely general Higgs potential, perform a GL{N,C) 
transformation in the space of doublets that brings the quadratic term to the form MoVo and 



then proceeds as discussed in Section |2.2[ Of course, the potential will contain very many 
different terms, so that the orbit space becomes a highly non-trivial multi-dimensional shape. 
However comprehending it is only a human limitation, and a hypothetical computer algorithm 
could be able to analyze this shape looking for edges, cusps and vertices. 

This method becomes much more useful when the number of distinct terms becomes small. 
In particular, when the dimension of the neutral orbit space is three or less, the shape can be 
relatively easily visualized, and one can develop a much more intuitive picture of the model 
than from the usual algebra. For example, in the 3HDM, this situation takes place for the 
following finite Higgs-family symmetry groups (based on the classification of [I3j): A4, S4, 
A(54)/Z3, and S(36) (the last two cases not discussed in this paper). It would be interesting 
to see if other useful examples appear in other models. 



4.5 Conclusions 

In summary, we have presented a simple yet powerful and very intuitive geometric approach to 
minimization of highly symmetric potentials in non-minimal Higgs mechanisms. This method 
is capable of giving the positions of the global minima with very little calculations; in partic- 
ular, it avoids the need to differentiate the potential, solve for stationary points, and check 
the positivity of the hessian. In a single picture, it shows all points which can be the global 
minimum for any values of the parameters of the potential. 

For illustration purposes, we have applied this method to A^- and S'4-symmetric three- 
Higgs-doublet models and found all points of global minimum. By doing this, we have also 
proved that the vacuum point used in at least three recent phenomenological analyses of the 
A4 case was not the global, but only the local minimum. We have also observed an unexpected 
approximate symmetry linking Higgs spectra at different minima and discussed its origin. We 
believe that this method can become a useful tool in all situations where minimization of 
highly symmetric functions is required. 
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A Higgs spectra of the 6'4-symmetric potential 

In the case of 5'4-symmetric 3HDM we analyzed the two simple vev alignments in the main text. 
They were shown to be approximately related to each other by the unexpected symmetry of 
orbit space. The remaining two points also follow this pattern. The alignment (1, i, 0) becomes 
the global minimum if 

A2<0, |A2|>|A3|, Ai>A3, 4A0 + A3 > 3IA2I , (28) 
and at this point we have 



V 



4Ao + A3 - 3|A2 
and the neutral Higgs masses are 



, m^± = -|A2|^;^ -(|A2|-A3y, (29) 



mjj. — ^(Ai— A3)v^ (double degenerate) 

^(|A2|+A3)^^ ^(|A2|+A0^^ 
ml = ^(4Ao + A3-3|A2|)^;2 = -^Mo. (30) 

The alignment (±1, e*'^/^, e"*'^/^) is the global minimum if 

A2<0, |A2|>|Ai|, A3>Ai, 4A0 + A1 > 3IA2I , (31) 
and at this point we have 

2 A^/3Mo 2 l|Ai2 1/iAi A^2 rQo\ 

' = 4Ao + A,-3|A2| ' "^-t = 2 l^^l" ' i^l^^l - ' ^'2) 

with the neutral Higgs masses 

mjj. — (a + 6 ± VaP^+~^) (double degenerate), a — ^"^^^ , b — 



6 



4A0 + A1-3IA2I 2 2 



ml = g ^^-^ v^ = ^Mo. (33) 

Again, we observe the perfect Ai -H- A3 symmetry in v"^, in the minimum conditions and in 
the charged Higgs masses. 



B Higgs spectra of the A4-symmetric potential 

Let us write for completeness the Higgs mass spectrum at all four possible points of global 
minimum found in the main text. 
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The vev alignment (1, 1, 1) remains stable in the presence of non-zero A4 if it is not too 
large: 

A2<12A?, A2<2(A3 + |Ai|)(A2 + |Ai|). (34) 
The value of is the same as in (20), while the masses become 
1, 



m ± 



\K^\±—K,]v', m^ = -(Ao-|Ai|y 



12 



5|Ai| + 3A2 + 2A3 ± V(|Ai| + 3A2 - 2A3)2 + 12A 



(double degenerate) . 



Note that the presence of A4 splits the charged Higgs masses while it preserves the 
degeneracy of the neutral Higgses. 

The vev alignment (1, 0, 0) is also stable if A4 satisfies 

A2<4(Ai + |A3|)(A2 + |A3|). (35) 

The value of and the masses of the charged Higgs and the non-degenerate neutral 
bosons are the same as in (23), while the neutral Higgses from the second and third 
doublets get masses 
,2 



V 



Ai + A2 + 2IA3I ± ^ - ^2? + A 



(double degenerate) . 



Note that this Higgs spectrum remains 2HDM-like as it was in the 5*4 case. 

The alignment (1, e*", 0) for a generic a parametrizes the points around the rim of the 
cone. For a given value of A4, the value of a corresponding to the global minimum is 
fixed by the relation 

A4 Ai - A2 , , 

sin 2a = , cos 2a = , (36) 

v/(Ai - A2)2 + AI ' V(Ai - A2)2 + AI 

Geometrically, this result means that the global minimum lies on the rim in the same 
direction as the vector n projected on the plane of the rim. In this case we obtain 

2 4y3Mo 



where 

A = -(Aic^ + A2S^ + A4C„s, 



4Ao + A3 - 3A 
1 



(Ai-A2)2 + Al-(Ai + A2) 



> 



(37) 



(38) 



and the mass spectrum is 

1 



v'^k, -t;'(A-A3) 

2 ' 4 V AI 



m 



ml 



-(A3 + A) + (1 ± cos 3a) •y/(Ai - Aa)^ + Aj 



-t;2(A, + A2 + 2A) 



(39) 
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Note that presence of cos 3a in one pair of masses is natural and it reflects the triangle 
symmetry of the A4 orbit space shown in Fig. |3] All other quantities are rotationally 
invariant, corresponding to the rotational symmetry of the cone. 

In the limit A4 — > 0, we get minimum at a — )■ tt/2, A — > — A2, and these spectra turn 



into the S'4-spectra found in Section 3.4 Note also that the three points on the rim with 
cos 3a = ±1 can never be "good minima" because we get mass terms with coefficients 
±(A3 + A), which cannot be made positive simultaneously. These three points lie, in 
fact, on the three long directrices of the cone. 

Finally, in the case of the alignment (±1, e*", e~*"), with a = 7r/3, we introduce conve- 
nient notation 



Aa 



sm7 



v/(Ai-A2)2 + A|' 



cos 7 



Ai - A2 



v/(Ai-A2)2 + A2 



and 



A = Aic2 + A^sl + A^Co^Sa 



Then, the value of v'^ is 



the charged Higgs masses are 



Ai + A2 + cos(2a - 7)\/(Ai - Aa)^ + A; 



An + A 



(40) 



(41) 



(42) 



and - -v"^ Ai + A2 + -^A, 
4 V 73 



(43) 



The neutral Higgs spectrum contains, as usual, h with mass = 2Mo/a/3 and two 
pairs of degenerate Higgses with masses 



m 



Hi 



V 

6" 



A3 + -(A1 + A2) -4A 



± ^/ ( A3 - -(Ai + A2) + 2A ) +3 [(Ai - A2)2 + AI] sin2(2a - 7) 



.(44) 



In the S'4-symmetric case (A4 = and 7 = 0), we recover the results of Section 3.4 



For the first three cases, the Higgs spectra were also explicitly written in [7]. We have checked 
that our results fully coincide with theirs. For the last case, the minimum 
was not explicitly written and analyzed by the authors of [7]. Instead they performed a 
numerical analysis of various minima of the type (r, e*°, e~*°), but their expressions, taken 
literally, become indeterminate at r = ±1 and a = 7r/3. It is conceivable that this last point, 
if treated as a special case, can still be recovered from their starting expressions. 
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